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Abstract. - Using exact diagonalization, Monte-Carlo, and mean-field techniques, characteristic 
temperature scales for ferromagnetic order are discussed for the Ising and the classical anisotropic 
Heisenberg model on finite lattices in one and two dimensions. The interplay between nearest- 
neighbor exchange, anisotropy and the presence of surfaces leads, as a function of temperature, to a 
complex behavior of the distance-dependent spin-spin correlation function, which is very different 
from what is commonly expected. A finite experimental observation time is considered in addition, 
which is simulated within the Monte-Carlo approach by an incomplete statistical average. We find 
strong surface effects for small nanoparticles, which cannot be explained within a simple Landau 
or mean-field concept and which give rise to characteristic trends of the spin-correlation function 
in different temperature regimes. Unambiguous definitions of crossover temperatures for finite 
systems and an effective method to estimate the critical temperature of corresponding infinite 
systems are given. 



Introduction. — The theory of collective magnetic 
order is usually concerned with infinite systems [iHl] . Real 
magnetic samples, however, have a finite size, and mag- 
netic properties are measured during a finite observation 
period. This gives rise to interesting questions as a mat- 
ter of principle. For magnetic nanoparticles it also has 
profound physical consequences. The recent advances in 
controlling and measuring magnetic properties of nanopar- 
ticles [S] as well as applications for magnetic data storage 
technology [6] rely on the fact that information, i.e. the 
magnetic state of a finite, small area representing a single 
bit, is stable over a finite observation time. 

The Curie temperature of a ferromagnetic sample is a 
typical concept that must be reconsidered for nanosizcd 
objects. From the experimental point of view, this is a 
well-defined quantity which can be measured, e.g. as a 
function of the size L of the nanoparticle [7l|8| . The mag- 
netic susceptibility and the specific heat stay finite but 
show enhancements at T = Tc{L), which defines a "Curie 
temperature" up to some residual arbitrariness. 

From the theoretical point of view, there is no Curie 
temperature as there is no phase transition, and actually 



not even a concept of a thermodynamic phase, in a system 
with a finite number of degrees of freedom. Nevertheless, 
one would like to define Tq (L) roughly to be the tempera- 
ture where the ferromagnetic alignment of the spins within 
the particle becomes stable against thermal fluctuations. 
In a simple mean- field picture, one finds Tq[L) < Tc(oo) 
where Tc(oo) is the precisely defined Curie temperature 
of the corresponding (infinite) bulk system. 

The transition, or smooth crossover, from a param- 
agnetic (PM) to an ordered state in the nanosystem 
must be described within models of interacting micro- 
scopic (atomic) spins [HITO]. Below Tc{L) the different 
microspins are tightly bound together and form a huge 
macrospin [51110] . For a finite-size and spin- isotropic sys- 
tem, the direction of the macrospin fiuctuates strongly, i.e. 
the magnetic state of the system is not stable temporally. 
Anisotropics give rise to superparamagnetic (SPM) behav- 
ior for temperatures above the so-called blocking temper- 
ature Tb = Tb(L), i.e. for Ty,[L) < T < Tc{L). Both 
quantities are size-dependent. It is rather the blocking 
temperature Th{L) than Tc{L) that is relevant for stor- 
age technology since it characterizes the crossover from 



p-1 



F. Author et al. 



the stable ferromagnetic (FM) state at low T to the SPM 
state where the system "switches" between different en- 
ergy minima determined by magnetic anisotropies. 

The topic is complicated by the fact that Tb(-^) cannot 
be considered as a pure property of the system. It must be 
seen as a relative value, which depends on the observation 
time r. For r — ^ oo, there is no blocking of the mag- 
netization as eventually the anisotropy energy barrier is 
overcome by thermal fluctuations or even due to quantum 
tunneling, and hence T\^{L) — >■ 0. For r — > (referring 
to e.g. laser-probe methods), T]^{L) Tc{L) while for 
intermediate r (like in spin-polarized scanning tunneling 
microscopy [H]) < rb(L) < Tc(L). 

For the infinite system, a magnetic phase transition is 
characterized by a divergence of the correlation length ^, 
which characterizes the spatial decay of the spin-spin cor- 
relation function. For a nanosizcd system, on the other 
hand, it is by no means clear how the crossover at the 
Curie temperature Tq (L) manifests itself in the correlation 
function. The purpose of the present paper is to provide a 
systematic study of the spin-spin correlation function for 
isotropic and anisotropic classical spin models in different 
dimensions with a finite and, in a thermodynamic mean- 
ing, small number of microspins. To this end, we compare 
results obtained by Landau mean-field theory, exact diag- 
onalization, and Monte-Carlo data. As our main result we 
find that three different temperature scales, the blocking 
temperature T\^{L) and the Curie temperature Tc{L) of 
the finite spin system as well as the Curie temperature of 
the infinite bulk rc(oo), can be read off from a suitably de- 
fined spin correlation function, which is accessible to scat- 
tering experiments. We propose a simple three-parameter 
fit formula for the correlation function, which turns out to 
be very effective in describing the numerical data for the 
entire temperature range and may serve to give a definition 
for the Curie temperature that is consistent with the usual 
estimates of Tc{L) based on the magnetic susceptibility or 
the specific heat. Finally, the blocking temperature, rel- 
ative to the observation time, can easily be accessed by 
interpreting the Monte-Carlo sweeps as time steps. Our 
results demonstrate that studies based on microscopic spin 
models in the superparamagnetic regime must go beyond 
the mean-field level. 

Spin correlation function. The correlation func- 
tion between two spins Si and Sj at sites and Tj is 
given by 

G(r) = (5,5,), (1) 

where (• • •) is the canonical thermal average at tempera- 
ture T. For a translationally invariant bulk system, the 
correlation function is homogeneous and depends on the 
translation vector r = — r only, while in case of a finite 
system it depends on r and on the reference site in addi- 
tion. For the following discussion we define an averaged 
correlation function which is independent of the direction 



and depends on the distance r = |r| only: 

\ri—rj \ —r 

Here, the sum in the first term runs over all n(r) pairs 
separated by the distance r. G{r) directly refers to X- 
ray- or neutron-scattering experiments. Furthermore, we 
define a "connected" correlation function 

G(r) = G(r) - , (3) 

where AI = \ ^^"^ ^ being the number of sites. 

Apart from a constant factor. Si is the local magnetic mo- 
ment at site ri, and thus M is the magnetization of the 
nanosystem. For temperatures above the blocking tem- 
perature, M averages to zero. But even below Tb{L) the 
magnetization vanishes, M = 0, in an exact calculation. 
The reason for the vanishing magnetization is the infinitely 
long "observation time" or time averaging in exact calcu- 
lations. If, on the other hand, the average (• • •) is inter- 
preted as a time average and if the time constant defining 
the average (life-time of a state) is large as compared to 
the observation time r, we have a finite M for tempera- 
tures below the blocking temperature corresponding to r. 
Hence, the function G(r) will strongly differ from G(r) or 
even vanish. 

Ising model. — We start the discussion by consider- 
ing the ferromagnetic (J > 0) Ising model: 

H^-JY,S^zSJ,, (4) 

on a one-dimensional chain of length L with open bound- 
ary conditions. The sum runs over all pairs of nearest- 
neighbor sites. The average correlation function G(r) is 
easily calculated analytically. Results for L = 10 and dif- 
ferent temperatures are displayed in Fig. ([TJi). We note 
that there is a simple exponential decay of the correla- 
tions. G{r) (X exp(— r/^), on a length scale ^ which at 
low temperatures exceeds the system size. Although in 
the low-temperature regime the microspins are perfectly 
aligned ferromagnetically and although they become un- 
corrclated on length scales much smaller than the system 
size in the high-temperature limit, there is no meaningful 
"Curie point" that could be extracted from G(r). Qual- 
itatively, not much happens as a function of T. The ab- 
sence of a Curie point, Tq{L) = 0, is of course not unex- 
pected. It corresponds to a featureless magnetic suscepti- 
bility x{T) and to the fact that Tc(oo) = for the infinite 
Ising chain. Only at T = 0, the system "freezes" in one 
of the two ferromagnetic ground states, i.e. M = 1, and 
therefore, due to our definition ([3]), the correlation func- 
tion discontinuously jumps to G(r) = 0. This might be 
expressed as a vanishing blocking temperature, T\j{L) ~ 0, 
reflecting the fact that an exact calculation corresponds to 
an infinite observation time t. 
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Fig. 1: Spin correlation function G(r) (see Eq. ((3]), symbols) 
calculated for an open Ising chain consisting of L = 10 sites 
(a) and for an open 5x5 Ising square lattice (b). The data 
are fitted with Eq. ([5]). The temperatures Ti (a) and Ti, T2, 
Is (b) are highlighted by thick red lines. 



For a finite two-dimensional Ising array with L = 5 x 5 
sites, the situation changes completely. Calculations for 
L = 25 are easily done by numerically exact diagonaliza- 
tion. Results obtained for the average correlation function 
G(r) [Eq. ([3])] are shown in Fig. [T}d. The result is surpris- 
ing: We find two or, including T = (see below), three 
different crossover temperatures. 

For high temperatures, the correlations decay exponen- 
tially, see T = 2.8J, for example. Below a temperature 
Ts, however, the trend can no longer be fitted by an expo- 
nential of the form exp(— r/^). We find T3 w 2.2J. This is 
close to the bulk Curie temperature of the two-dimensional 
Ising model T3 w Tc(c») = 2/ ln(l -h v^) J « 2.27J. Upon 
lowering T we then find another temperature T2 ~ 1.9 J 
which is characterized by a change of the curvature of 
G{r). Below T2, the trend of G{r) is no longer convex but 
concave until r hits the system boundary. As the exact- 
diagonalization data correspond to an infinite observation 
time, the third temperature scale Ti is trivially given by 
the vanishing blocking temperature Ti ~ Tij{L) ~ 0. 

Model correlation function. To extract the dif- 
ferent temperature scales T3,T2,Ti, we propose to fit the 
correlation function to the following expression with three 
temperature-dependent parameters: 
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G(r) = G(r, T) « B{T)e-''''^^^ + y{T) . 



(5) 



Fig. 2: Spin correlation function G(r) (see Eq. ([T]), sym- 
bols) calculated for a 20 x 20 Ising square lattice with open 
boundaries starting from an edge (left) and from the cen- 
tral site (right) as obtained by Monte-Carlo simulations at 
T = 0.5 J ^ rc(oo). The lines are guides to the eyes only. 



In this way is defined by the temperature where y(T) 
becomes finite, i.e. where a deviation from a purely ex- 
ponential decay of G{r) is found. For the infinite system, 
L — > 00, this happens right at T3 = Tc(oo) where a power- 
law decay is expected [UIH]. Thereby, the simple fit for- 
mula, Eq. ([5]), will provide a rough estimate for Tc(oo) 
based on a single calculation of G(r) for a finite system. 
Our Monte-Carlo calculations show that this estimate re- 
liably gives Tc(oo) within an error of the order of 1%. For 
example, using the fit for Monte-Carlo data obtained for 
the 5x5 Ising array, we find = 2.20 ± 0.02 while for 
the 8x8 lattice = 2.26 ± 0.02. This represents a cheap 
but rough way to get Tc(oo) from a slow annealing of a 
single finite system. For temperatures below the block- 
ing temperature, we have G(r) = (G(r) = 1). Hence, 
Ti ^ T^{L) is indicated by B{T) = and y(T) = 1 when 
fitting the data using Eq. ([5]) . The main purpose of Eq. 
([5]), however, is to get an estimate for T2 where G(r) is 
a linear function to a good approximation. A linear G(r) 
requires e{T) — > 00 for T ^ T2. To get a finite slope, we 
also need B(r) -> 00 for T ^ T2, and finally y{T) 00 
since |G(r)| < 1. 

We have used Eq. (P to fit the unknowns £(T), B{T) 
and y{T) to numerically exact data for one-, two- and 
three-dimensional Ising systems of different size L as well 
as for finite isotropic and anisotropic Heiscnberg systems. 
It turns out that the quality of the fit is exceptionally good 
in the entire temperature range, sec the lines in Fig.[T^ and 
b, for example. For temperatures Tb{L) < T <T2 wc find 
the concave trend with negative values for B{T) and e{T) 
while y{T) > 0. For T2 < T < Tc(oo), both B{T) and 
e{T) are positive while y{T) < 0, and the trend of G(r) is 
convex. Finally for T > T3 both B{T) and e{T) remain 
positive while y{T) vanishes leaving an exponential decay 
of G(r). 

Curie temperature. — We identify T2 with the 
"Curie temperature" of the system, T2 ~ Tc{L). This 
provides a meaningful definition of Tq{L) for a finite sys- 
tem that is based on the spin correlation function. It is 
motivated by the physical idea that at Tc(L) the correla- 
tion length exceeds the system size but additionally takes 
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into account that the nanosystem is bounded by surfaces. 
For an infinite system the spin correlation function is al- 
ways convex, i.e. its slope is negative but increasing as 
a function of increasing distance r. The unusual concave 
trend of the correlation function at lower temperatures 
must therefore be a direct consequence of the presence of 
surfaces. This is demonstrated with Fig. [2] for a larger 
system consisting of 20 x 20 Ising spins on a square ar- 
ray. Due to missing nearest neighbors at the nanoparticle 
surface, fluctuations of the local spins arc stronger and re- 
sult in a reduced average surface magnetization. This also 
implies a strongly decreasing correlation function G{r) 
(see Eq. ([ij) close to the surface and along any direction. 
For the averaged correlation function G(r) (see Eq. ([2])), 
this surface effect competes with the convex bulk trend 
of G(r). The surface effect dominates for T > Tq{L) and 
drives the nanosystem to a "paramagnetic" state while for 
T < Tc{L) the bulk of the nanosystem causes an ordered 
superparamagnetic state, and the surface manifests itself 
in stronger fluctuations of the spins and a concave trend 
of G(r) only. Note that the fit with the model correla- 
tion function, Eq. ([5]), allows to characterize Tc{L) by a 
divergence of the parameter £(T), which therefore might 
be called a "virtual correlation length". On the other 
hand, the Curie temperature of the infinite system Tc(oo) 
is given by a divergence of = —e{T)hi[—y{T)/B{T)] 
since y{T) becomes finite at rc(oo) in the fit. is the 
distance at which G{r) vanishes, G{r^) = 0. The distance 
is "virtual" because it is always larger than the system 
size (see Fig. 

Blocking temperature. — The question of a finite 
blocking temperature Tb can be addressed by Monte- 
Carlo simulations when interpreting Monte-Carlo sweeps 
as time steps [13]. A finite number of Monte-Carlo 
sweeps corresponds to an incomplete statistical average 
and thus to a finite observation time r. Calculations have 
been performed for the classical Heisenberg model with 
nearest-neighbor exchange J > and uniaxial magnetic 
anisotropy K > Q on finite two-dimensional square lat- 
tices: 



h = -jY^s,s,-kY,{s;] 



(6) 



Each Monte-Carlo run consists of up to 200 temperatures 
with up to lO'' sweeps per temperature. 

An example for K = 0.6 J is given in Fig. [3J:. The fit of 
the Monte-Carlo results by Eq. ^ is accurate for all tem- 
peratures and all values of K such that the three different 
temperature scales, T\^(L),Tc{L), and Tc(oo) can be ex- 
tracted easily. The finite blocking temperature manifests 
itself in the jump of G(r) as a function of T which is due 
to the jump of M at T « 0.6J. The order of magnitude 
for Tb seems to be given by K. However, T\j{L) decreases 
with increasing observation time r, i.e. with an increasing 
number of Monte-Carlo sweeps per temperature. At fixed 
L, we find T]^{L) — !• logarithmically if r — )■ oo. 
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Fig. 3: Model G(r) [see Eq. ([5}] as obtained by fitting the pa- 
rameters to results of Monto-Carlo simulations (10® sweeps per 
temperature) for an anisotropic Heisenberg model {K — 0.6J) 
on an 8 X 8 square lattice with open boundary conditions, (a) 
Temperature dependence of the virtual distance r„ (see text), 
(b) Virtual correlation length e(T). (c) G{r) as a function of r 
and T. The temperatures cxj > Tc(oo) > Tc{L) > Ty,{L) > 
are highlighted. In each of the four corresponding temper- 
ature ranges, the typical Monte-Carlo time dependence of 
M = M{t) (X C^-Si) is shown in insets (see text for dis- 
cussion). 



Fig. [3J: nicely demonstrates that an astonishingly com- 
plex behavior of the spin correlation function is found for 
finite anisotropic nanosystems. The qualitatively differ- 
ent physics within the different temperature ranges, i.e. 
< 'rb(L) < Tc(i) < Tc(oo) < oo, also shows up in the 
qualitatively different behavior of the order parameter M 
as a function of (Monte-Carlo) time during the simulation, 
see insets in Fig. 131:;: Below T\y{L), on the scale set by the 
observation time, the magnetization freezes in one of the 
values corresponding to the degenerate energy minima of 
the anisotropic model; for T'b(L) < T < Tc{L) the mag- 
netization switches between these values with a switching 
time which is much smaller than the "magnetic lifetime" of 
a state; for Tc{L) < T < Tc(oo) the system still switches 
but the lifetimes are comparable to the switching times; 
and finally above rc(oo) correlations decay exponentially 
and M ^ 0. 
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Fig. 4: Curie temperature, as obtained from tlie maximum 
of tlie specific fieat (open squares), from tiie maximum of tlie 
susceptibility (open circles) and from the correlation-function 
fit (filled circles), as functions of the (inverse) linear system 
size, and bulk Curie temperature (filled squares), as obtained 
from the fit of the correlation-function at the respective sys- 
tem size. Solid line: (Tc{oo) - Tc{L))/Tc{oo) = [L/Lo)-^''^ 
corresponding to the exponent v = 1 and with L, 
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Dependence on the system size. Note that this 
physics is characteristic of a finite system: For constant r 
but increasing system size L — >■ cxj all three temperatures 
merge, and we are left with a single critical temperature 
only, the Curie temperature of the infinite system Tc(oo). 
With increasing L, but keeping the temperature fixed, the 
curvature of G(r) increases, i.e. it becomes less convex, 
changes from convex to concave, or becomes more con- 
cave. This is due to the less and less important effect of 
the nanoparticle's surfaces. At the same time G{r) and 
its slope increase. This implies that Tq{L) is an increas- 
ing function of the system size. The same holds for the 
blocking temperature since with increasing L the energy 
of the anisotropy barrier increases and higher temperature 
is needed to induce a thermal switching of the magnetiza- 
tion. 

The system-size dependence of the Curie temperature, 
as obtained from the fit of the correlation function, is dis- 
played in Fig.(|4]). Comparing Tc{L) = T2 with the Curie 
temperature defined by the maximum of the magnetic sus- 
ceptibility and by the maximum of the specific heat, we 
find that the asymptotic behavior is approached signifi- 
cantly faster when using the correlation- function fit. The 
latter also provides a reliable estimate for the bulk Curie 
temperature as is also shown in Fig. (|4]). In addition. Fig. 
([5|) shows the dependence of the blocking temperature on 
the system size. For fixed system size and with increasing 
Monte-Carlo time, i.e. with increasing number of sweeps, 
the blocking temperature decreases logarithmically (see 
data for L = 8 X 8 in Fig. (P). 
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Fig. 5: Blocking temperature, as obtained from the 
correlation-function fit, as a function of the (inverse) linear 
system size for the D — 2 anisotropic Heisenberg model with 
K = 0.6J and using 1.0 • 10" Monte-Carlo sweeps. Lines con- 
nect the data points. For the L = 8 x 8 system, the dependence 
of Tb on the Monte-Carlo time is given by the red circles as in- 
dicated. 



Our analysis shows that Tq{L) satisfies the finite-size 



scaling law (Tc(oo) -rc(L))/Tc(oo) = (L/Lq 



-1/-DA, 



[T5] . v'Lo corresponds to a microscopic length scale, and 
its order of magnitude is one. The shift exponent A^. is 
related to the exponent of the correlation function via = 
\/v. In case of the anisotropic Heisenberg model, our 
data for Tc{L) for system sizes up to L = 19 x 19 = 361 
are consistent with ly = 1.0 (and ^/Lq = 1.15). This is 
different from the classical exponent {i/ ^ 0.5) and agrees 
with the exponent for the D = 2 Ising model [v = 1). 

Note, that our results for the size dependence of the 
blocking temperature are also consistent with the same 
scaling law that describes Tc(L), see Fig. ([S]). The ex- 
ponent As, however, is different and slightly larger than 
1/v for V = \. Actually, it is by no means clear that the 
blocking temperature should satisfy a scaling law since for 
the infinite system it has no meaning independent from 
the Curie temperature. This point is beyond the scope of 
the present study but deserves further investigations. 

Static mean-field theory. The characteristic 
trends of the spin-correlation function in the different tem- 
perature regimes are strongly determined by the presence 
of surfaces, see Fig.[2j It is tempting to simply explain the 
concave trend of the correlation function at temperatures 
below Tc{L) by the reduced coordination numbers at the 
nanoparticle surfaces. We have checked this by performing 
calculations using periodic boundary conditions. In fact, 
a convex curvature of the correlation function is found in 
one and in two dimensions and for all temperatures in 
agreement with previous work [3]. Since the surface-to- 
volume ratio is smaller in one as compared to two dimen- 
sions, the simple coordination-number argument is also 



p-5 



F. Author et al. 



consistent with the absence of a Curie temperature in one 
dimension. We have checked our findings by performing 
corresponding calculations for three-dimensional finite lat- 
tices: In fact, the low-tcmpcraturc concave trend is found 
to be even more pronounced here. 

Simple coordination-number arguments arc included in 
the Landau theory of magnetic systems bounded by sur- 
faces [T6HI9] where the usual Landau free-energy func- 
tional in considered but with an additional surface free- 
energy term. The resulting Landau mean-field theory is 
essentially equivalent to static mean-field theory for a dis- 
crete spin model. For a finite system, this is easily im- 
plemented numerically, and we have evaluated G{r) for 
different one- and two-dimensional lattices studied here. 

As expected, static mean-field theory gives a phase tran- 
sition rather than a smooth crossover. The mean-field 
Curie temperature J'^'^^' of the finite system very much 
depends on the size and the geometry of the underly- 
ing lattice. If interpreted as Tc{L), the mean-field Curie 
temperature T^^'^ yields a strong overestimation. Here, 
however, our question is whether besides T^^'^ there is a 
crossover temperature at which the spin-correlation func- 
tion changes qualitatively from convex to concave. 

The spin-correlation function can be obtained in two 
ways, either directly by computation of the thermal av- 
erage {SiSj) or as the response of the local magnetic 
moment at site i to a local magnetic field at site j, i.e. 
d{Si) /dBj. In principle, both ways are equivalent because 
of the fluctuation-dissipation theorem (SiSj) — (Si) (Sj) = 
Td{Si) /dBj. This is not respected by static mean-field 
theory, which just neglects non-local correlations. There- 
fore, spin correlations can be addressed via the linear- 
response relation G(r) = Td{Si) / dBj only. For an infinite 
translationally invariant lattice, this yields the Ornstein- 
Zernike form for the correlation function. Here, for fi- 
nite systems, we determine (Si) numerically by solving 
the static mean-field equations and compute the derivative 
with respect to Bj numerically. Looking at the resulting 
averaged correlation function G'(r), we always find a con- 
vex trend, for any system size and dimension. This shows 
that our results and the crossover temperature Tc{L) can- 
not be captured by a mean-field or Landau approach and 
therefore represent a correlation effect for which simple 
coordination-number arguments must be taken with care. 

Conclusion. — As compared to infinite bulk systems, 
the theoretical description of collective magnetic order is 
more involved for nanosized materials. Due to the finite 
system size there are no clear-cut regions in parameter 
space where ferromagnetic order is realized. Furthermore, 
the magnetic state is not stable temporally and conse- 
quently the order parameter, i.e. the magnetization of the 
nanoparticlc, fluctuates with a time constant that has to 
be compared with the (experimental) observation time. 
These facts give rise to ambiguities in the deflnition of the 
Curie temperature and imply the existence of a second 
temperature scale, the blocking temperature, which again 



cannot be deflned precisely. 

Our studies based on different analytical and numerical 
techniques have demonstrated, that a meaningful defini- 
tion of the Curie temperature of a finite spin system can 
be given that relies on the analysis of a suitably defined 
average spin correlation function. Upon lowering the tem- 
perature, the correlation function changes its curvature at 
Tc(L). This definition is consistent with the usual con- 
cepts and comes closest to our expectation that the Curie 
point is the temperature at which the correlation length 
exceeds the particle size. In particular, it accounts for the 
delicate interplay between the (bulk) tendency to order- 
ing and the (surface) tendency to enhance fluctuations. 
We have shown that the concept can be applied to differ- 
ent one-, two- and three-dimensional classical spin models 
and that Tc{L) can be extracted with an accuracy that 
even allows to determine a shift exponent. 

We could further demonstrate that the blocking tem- 
perature scale is accessible with a Monte-Carlo approach 
by performing an incomplete statistical average. A sharp 
jump is visible in the average spin-correlation function at 
Tb(i). However, the blocking temperature is defined with 
respect to an observation time only (a finite number of 
Monte-Carlo sweeps) and therefore represents a relative 
quantity. 

Concluding, the combined application of exact diago- 
nalization, Monte-Carlo, and mean-field techniques has 
uncovered a strikingly complex behavior of the spin cor- 
relations in nanoparticlcs with qualitatively different tem- 
peratures scales. Since simple coordination-number argu- 
ments are unable to give a quantitively correct picture of 
the physics, the temperature trends must be seen as effects 
of strong spin correlations. 

* * * 

Support by the DFG (SFB 668, projects B3, B4, A3, 
and A14) and by the Cluster of Excellence "Nanospin- 
tronics" Hamburg is gratefully acknowledged. 

REFERENCES 

[1] Binder K., Rep. Progr. Phys., 60 (1997) 487. 
[2] Fisher M. E., Rev. Mod. Phys., 46 (1974) 597. 
[3] Fisher M. E., Barber M. N., Phys. Rev. Lett, 28 (1972) 
1516. 

[4] Wu T. T., McCoy B. M., Tracy C. A., Barouch E., 

Phys. Rev. B, 13 (1976) 316. 
[5] Krause S., Herzog C, Stapelfeldt T., Berbil- 

Bautista L., Bode M., Vedmedenko E.Y., and 

WiESENDANGER R., Phys. Rev. Lett, 103 (2009) 127202. 
[6] Service R., Science, 314 (2006) 1868. 
[7] Nepijko S. A. and Wiesendanger R., Semiconductor 

Physics, Electronics and Optoelectronics, 2 (1999) 5. 
[8] KOTZLER J., GORLITZ D., KURFISS M., VON Sawilski L. 

and Vedmedenko E. Y.„ Phys. Rev. B, 73 (2006) 224425. 
[9] Landau D. P., Phys. Rev. B, 13 (1976) 2997. 



p-6 



Title 



[10] DiMlTROV D. A. and Wysin G. M.„ Phys. Rev. B, 54 
(1996) 9237. 

[11] WiESENDANGER R., Rev. Mod. Phys., 81 (2009) 1495. 
[12] Campanino M. and Ioffe D. and Velenik Y.„ Euro- 

phys. Lett, 62 (2003) 182. 
[13] E. Y. Vedmedenko, Phys. Stat. Sol. (b), 244 (2007) 

1133. 

[14] M. N. Barber, m: Phase Transitions and Critical Phe- 
nomena, edited by C. Domb and J. L. Lebowitz, (Academic 
Press, New York, 1983) 

[15] O. IGLESIAS and A. Labarta, Phys. Rev. B, 63 (2001) 
184416. 

[16] Mills D. L., Phys. Rev. B, 3 (1971) 3887. 
[17] Binder K. and Hohenberg P. C, Phys. Rev. B, 6 (1972) 
3461. 

[18] Kumar P., Phys. Rev. B, 10 (1974) 2928. 
[19] LUBENSKY T. C. and Rubin M. H., Phys. Rev. B, 11 
(1975) 4933. 



p-7 



